Towards Heisenberg Limit in Magnetometry with Parametric Down Converted 

Photons 
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Recent theoretical and experimental papers have shown how one can achieve Heisenberg limited 
measurements by using entangled photons. Here we show how the photons in non-collinear down 
conversion process can be used for improving the sensitivity of magneto-optical rotation by a factor 
of four which takes us towards the Heisenberg limit. Our results apply to sources with arbitrary 
pumping. We also present several generalizations of earlier results for the collinear geometry. The 
sensitivity depends on whether the two-photon or four-photon coincidence detection is used. 



PACS numbers: 00000 



I. INTRODUCTION 



Parametric down conversion is a process that is used to 
produce light possessing strong quantum features. Pho- 
ton pairs generated by this process show entanglement 
with respect to different physical attributes such as time 
of arrival Q and states of polarization 0] . They are in- 
creasingly being utilized for very basic experiments to 
test the foundation of quantum mechanics and to do 
quantum information processing 0. S E| ■ It is also rec- 
ognized that entangled photon pairs could be useful in 
many practical applications in precision metrology in- 
volving e.g. interferom etry 0,0 013) imaging [a Hof. 
lithography [TllfT^. Il3l Hif and spectroscopy . There 
is a proposal [16| to use electromagnetic fields in NOON 
states to improve the sensitivity of measurements by a 
factor of N. Some implementations of this proposal exist 
|17|. In particular, the use of photon pairs in interfer- 
ometers allows phases to be measured to the precision 
in the Heisenberg limit where uncertainty scales as l/N 
|l8j as compared to the shot noise limit where it scales as 
I/Vn. This means that for large number of particles, a 
dramatic improvement in measurement resolution should 
be possible. 

In this paper we present an analysis of how parametric 
down converted photons could be very useful in getting 
better spectroscopic information about the medium. We 
demonstrate how the improvement in magneto-optical ro- 
tation (MOR) of light could be realized by employing two 
different schemes with collinear and non-collinear down 
conversion geometry in compared to use of coherent light. 
We calculate the resolution that can be achieved in the 
MOR's both by use of coherent light and down converted 
light. We discuss that the Heisenberg limit 0] could 
be reached in magnetometry by the use down converted 
light. 



II. MOR USING COHERENT LIGHT SOURCE 

Consider a single mode coherent light travelling in 
the z-direction and a linear isotropic medium made 



anisotropic by the application of the magnetic field B 
in the z-direction. The incident field can be written in 
the form 

E(z, t) — exp(— iujt + ikz)(xe x + y£ y ) + c.c. (1) 

The medium is described by the frequency and magnetic 
field dependent susceptibilities x± (^-0 ■ That means hori- 
zontally and vertically polarized components of the inci- 
dent light will rotate on travelling the medium of length 
I and the field at the exit can be written as 

E(7, t) — exp(— tut + ikl){xe x i + ye y i) + c.c. (2) 

The rotation of the horizontal and vertical components 
can be expressed by the relations 



Syl 



= R 



where 
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R = e w +e l * 

6 = M(x+-X-), 
0+= Mx+- 



sin ■ 



(3) 



(4) 

(5) 
(6) 



The corresponding quantum mechanical description 
can be obtained by replacing the classical amplitudes e x 
and e v by the annihilation operators a x and a y respec- 
tively. For measurements with coherent sources one can 
look at the intensities of the x and y components of the 
output when the input is x polarized with coherent state 
amplitude a x (See Fig. ^a). Then the measured quan- 
tities will be 



Ixl = (a x i f a x i) = \a x \ 2 cos 2 -, 



yl 



{a v i ] a y i) 



|a x | 2 sin 2 -. 



(7) 
(8) 



One can estimate the minimum detectable rotation an- 
gle 9 m by looking at the fluctuations AiV<j in the photon 
number difference between horizontal and vertical pho- 
tons, where the number difference operator is given as 



2 




FIG. 1: The setup for the Magneto-optical rotation of light 
by using (a) coherent source, type-II PDC photons with (b) 
collinear and (c) non-collinear geometry. 

Nd = dyV a yi ~ a x^o-xi- This expression is calculated to 
be (ATV^) 2 = | cv^ | 2 sin 2 # and since the fluctuation noise 
is 1 we obtain 9 m « 1/ yj (N) where (N) is the mean 
number of input photons which is equal to |a x | 2 . 

III. MOR USING COLLINEAR TYPE-II PDC 
AND TWO-PHOTON COINCIDENCE 

We now discuss how the results (|7|l and JSJ are mod- 
ified if we work with down-converted photons. We first 
consider the collinear case shown in Fig. C^b). The state 
produced in collinear PDC can be written by 

|Vw> = —t- V(-e^tanhr)» ff |n) y . (9) 
coshr 

n=0 

The value of the parameter r and the phase <p are re- 
lated to the pump amplitude of the nonlinear crystal that 
is used in the down conversion process and the coupling- 
constant between the electromagnetic field and the crys- 
tal. Note that the state \ipcoi) is a superposition of n pho- 
ton pairs of horizontally and vertically polarized modes. 
Inside the medium, these modes rotate with the same 
rotation matrix R given in Eq. Q): 

( am ) = R ( aH ). (10) 
\a V i J \ a v ) 

One can measure the intensity of each mode: 

Ih = (am^aHi) = sinh 2 r 

= (avi^avi) = ly. (11) 
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FIG. 2: The MOR plot of two-photon coincidence counts de- 
fined by the Eq. in collinear type-II PDC. r is the interac- 
tion parameter that defines the pumping strength used in the 
production of down converted photons and 9 = kl(x+ — X-)- 
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FIG. 3: The visibility of two-photon and four-photon counts 
defined by the Eqs. Il l 2ft and l|19|l respectively. 



And the two-photon coincidence count is: 

I H v = {aHV avv a H iavi) 

— cos 2 9 sinh 2 r cosh 2 r + sinh 4 r (12) 

Note the difference between Eqs. (JJJ and lfTS)l. With 
collinearly down-converted photons we measure a rota- 
tion angle that is twice as large compared with the angle 
for a coherent input. For r < 1 wc obtain the same 
result as given in |l5j . The fringe pattern and the vis- 
ibility is given in Figs. |21 and [3] One can calculate 
the minimum detectable rotation angle again by look- 
ing at the fluctuations in the photon number difference 
N d . This is given by (AiV d ) 2 = 4 sinh 2 r cosh 2 sin 2 9 
: I + (iV))(iV)sin 2 f5 « (7V) 2 sin 2 6» for large (N) where 
(N) = 2 sinh 2 r. Making (AN d ) ~ 1 we get 9 m m 
1/(N). Note that the sensitivity of this quantity is also 
improved by a factor of 1/ \/~N. 



3 




FIG. 4: (a) The normalized four-photon probability defined in 
Eq. 1171 , and (b) its envelope wrt the interaction parameter 
r in the non-collinear geometry. 



IV. MOR USING NON-COLLINEAR TYPE-II 
PDC AND FOUR-PHOTON COINCIDENCE 

Next, we discuss the non-collinear PDC case. We have 
found an arrangement shown in Fig. ^c) which is espe- 
cially attractive for improving sensitivity. The entangled 
photons are coming in two different spatial modes, a and 
b. While one mode (say a) is going parallel to B inside 
the medium, the other is going anti-parallel to it. At the 
exit we separate the H and V modes by polarizing beam 
splitters. The state of the input photons can be written 
in the form 20] 



IVv, 



cosh 



l 00 

-5- Y J V^+T(t&nhr) n \ip n ), (13) 



n=0 



where 



1 

(14) 

Here \m) av represents m vertically polarized photons in 
mode a. Inside the medium, " -I- " and " — " polariza- 
tion components of the modes a and b gain phases kl\+ 
and kl\- respectively. Thus we can write an effective 
Hamiltonian for the evolution of the state \ip non ) inside 
the medium as follows: 

H med = x+a+ ] a++X-a- a--x+b+ ] b + -x-b- ] b-, (15) 

where 

a ± = -^=(a H ±ia v ), b± = -^=(b H ± ib v )- (16) 
I 



The minus sign in front of the b± modes comes from the 
fact that they are travelling anti-parallel to the B field in- 
side the medium. Then one can calculate the probability 
of detecting four photons in each mode as: 



Pn 



= \{l aH l av lb H lb v \ exp(-itff me <i)|Vw)| S 



tanh r 9 . „. 

—rr- cos (20) 

cosh r 



(17) 



(b) 




FIG. 5: (a) The normalized four-photon probability defined in 
Eq. 1181 . and (b) its envelope wrt the interaction parameter 
r at the exit ports of PBS in the collinear geometry. 



where t is the duration for the state to evolve inside the 
medium. Note that this four-photon probability has the 
rotation angle that is four-times as large compared with 
the angle for a coherent input. The fringe pattern with 
respect to and the probability distribution with respect 
to r are shown in Fig. 01(a) and (b). 

Next we also examine the four-photon probability in 
the collinear case. The probability of finding two H- 
photons and two ^-photons at the exit ports of the po- 
larizing beam splitter is given by: 

Pcoi = \(2 aH 2 av \eiq>(-itH med )\ip col )\ 2 



tanh r 1 
cosh 2 r 16 



[1 + 3cos(260] 2 



(18) 



where = x+ a + a + + X- a - a - because of 

the collinear geometry. The normalized plot of this quan- 
tity with respect to the magneto-optical rotation angle 9 
and the envelope of the probability with respect to r are 
shown in Fig. [5] (a) and (b) . 

On the other hand one can also calculate the coinci- 
dence counts of four photons two-by-two at each detector 
as given by Glauber's higher order correlation functions: 



t I +2 t2 2 2\ 

j-hhvv = \o-hi avi am avi ) 

= (3 cos 2 6»- l) 2 sinh 4 rcosh 4 r + 4(3 cos 2 + 1) sinh 6 r cosh 2 r + 4sinh 8 r. (19) 

I 



The plot of this quantity for different values of the inter- 



action parameter r and the visibility are shown in Figs. 
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r = 0.1, scaled by 10 4 

r = 0.5, scaled by 5 

r = 1.3, scaled by 2x1 J 




FIG. 6: Four-photon coincidence counts defined in Eq. I|i9|l 
with different interaction parameter values in the collinear 
geometry. 



E]and|3| Note the distinction between Eqs. (jT5)l and (jl§|) 

which is a reflection of what the detector is set to mea- 
sure as we explain now. The former is the probability 
of the state |?/w(i)) to be projected onto the particular 
four-photon subspacc 1 22) i.e. Tr[|22)(22|p co /(i)] where 
Pcoi(t) = U\ip cc ,i)(if} co i\W and U is the unitary operator 
that represents the evolution of the state by the Hamilto- 
nian H rne( i in the collinear geometry. On the other hand, 
coincidence counting of four-photons at the detectors Dh 
and Dy (see Fig. ^b)) is represented by the expectation 
value (ajf t2 ov t2 «// 2 ay 2 ) = Tr [an^ay ^ 2 a.H 2 ay 2 p co i{t)]- 
Note here that the operator aH^ay^a^ay 2 has the 
spectral decomposition J^^m C nm \nm){nm\ and obvi- 
ously it contains the projectors of all (n + ?7j)-photon 
subspaces with nonzero coefficients G nm . Therefore four- 
photon counting process at detectors includes not only 
1 22) but all other states \nm) in \ip co i(t)). Here the 
state \nm) represents n and m photons in the an and 
ay modes respectively 



V. CONCLUSION 

We showed that the use of non-collinear type-II PDC 
light in MOR's increases the sensitivity by a factor of 



four in comparison to coherent light. We also give an 
argument that minimum rotation uncertainty scales to 
Heisenberg limit by the use of down converted photons. 
It should be noted that Heisenberg limit should be un- 
derstood as an approximate limit at a large mean photon 
number, that is, the rotation uncertainty approaches the 
order of 1/ (N) for large (N) [j| . The regime with an 
interaction parameter value of r = 1.3 has already been 
reached in the experiment pH giving entanglement of 12 
photons and an evidence also was given for entanglement 
up to 100 photons. 
APPENDIX A: FOUR-PHOTON PROBABILITY 

In this appendix we show the details of the calculation 
leading to the result given in Eq. (|T7|) . One can obtain 
the result first by solving the Schrodinger equation for the 
state \l aH loylfc^lfev) m the four-photon subspace of the 
electromagnetic field and having the inner product with 
the state \ip) n on- Since the parts of the Hamiltonian hav- 
ing a and b modes commute, we can solve the Schrodinger 
equation for the states \ l aH ^a v ) and \lb H ^b v ) separately. 
Let us start with a general time-dependent state in the 
an and ay modes which contain two photons totally; 

= c(t)|20) + d(t)\02) + /(t)|ll) (Al) 

with the initial condition |</>(0)) = |11). Solving the 
Schrodinger equation by using the effective Hamiltonian 
H = + X- a -^ a - gives us the result 



-itX\ 



1 



l V2 
cos(fit)|ll)] 



sin(Ot)|20) 



1 



V2 



sin(Q£)|02) 



(A2) 



where x — X+ + X- an d fi = X+ ~ X- ■ F° r a medium 
of length I, the angle fit corresponds to the MOR an- 
gle 9 which is given in Eq. The solution for the 
state \lb H lb v ) can be obtained just by replacing 9 by —9 
because the direction of propagation of the b modes are 
opposite to that of a modes inside the medium. This is 
the reason that the part of the effective Hamiltonian for 
the b± modes takes minus sign in Eq. (|15fl . Consequently 
we can write the solution of the Schrodinger equation for 
the state |l OH l av l bH l bv ) as: 



J 



exp(-itH medium )\l aH l av l bjI l bv ) = cxp(-iix) 

<g) exp(— itx) 

r 



-^sin(9|20) - -^sin<9|02) +cos0|ll) 
V 2 v 2 

— ^sin(9|20) + ^sin6»|02) +cos0|ll) . (A3) 
V 2 v 2 



Taking the inner product of this with the state \ip) n on in Eq. (|1 T|) . 
and having the absolute square gives us the result given 
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The result given in Eq. I|18|) can be obtained by fol- 
lowing the same method given above. 
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